﻿ Studiul semnalelor periodice 1 Scopul lucrării Se efectuează analiza în frecvenţă a semnalelor periodice şi se măsoară distorsiunile semnalelor de la ieşirea generatoarelor de semnale sinusoidale 2 Dezvoltarea în serie Fourier exponenţială a semnalelor periodice complexe Un semnal în timp continuu CR:→ x se numeste periodic daca exista o constanta strict pozitiva T astfel incat: (, Rxt T xt t+=∀∈ (1) )() Se numeste perioada a semnalului cea mai mica valoare a lui T care satisface egalitatea de mai sus ) Are loc urmatoarea proprietate: daca o functie ( tx este local integrabila si periodica de perioada T, atunci: TaT+ (), Rx t dt x t dt a=∀∈ (2) () 0 a∫∫ Deoarece integrala pe o perioada a functiei nu depinde de pozitia intervalului de integrare, de lungime egala cu perioada, pe axa timpului, aceasta marime se noteaza scurt cu: () dttx ∫ T Se spune ca un semnal periodic, de perioada T, satisface conditiile lui Dirichlet daca: tx este absolut integrabil pe o perioada: • () () ∞ ⋅=⋅=kccAaϕ 0 kkkkk {}() ,argsin2sin>⋅−=⋅−=kccAbϕ 0 kkkkk Cu referire la relatia (6): tx - marimea A se numeste componenta continua a semnalului ( 0) - semnalul () ω+⋅cos se numeste armonica de ordinul k a semnalului (tx, kktkAϕ0) 0>k kA k) - marimile se numesc valorile efective ale armonicilor (0> 2 - se numeste spectru de amplitudini reprezentarea numerelor A functie de k ,≥ 0 0kkω - se numeste spectru de faze reprezentarea marimilor ϕ kkω k functie de ,0≥ 0 4 Puterea semnalelor periodice ) Se considera un semnal oarecare (tx Daca exista marimea: τ 21 lim dttxP= () X∫ τ2 ∞→ τ τ − ) aceasta se numeste puterea semnalului ( tx ) Se poate arata ca, daca ( tx este un semnal periodic, de perioada T, atunci: 21(8) () dttxP= X∫ TT Au loc urmatoarele relatii Parseval: 2 22 ∞∞ (9) ⎛ 2⎞ kA +== ⎟⎜ ∑ 0∑ kXAcP −∞= kk 12= ⎠⎝ Prima egalitate are loc pentru orice semnal periodic, in timp ce a doua are loc doar pentru semnale reale Deoarece puterea unei sinusoide este egala cu patratul valorii sale efective, rezulta ca puterea unui semnal periodic real este suma dintre puterea componentei continue si puterile semnalelor sinusoidale care apar in dezvoltarea sa in serie Fourier armonica 5 Relatia intre coeficientii Fourier asociati unui semnal periodic si transformarea Fourier asociata restrictiei acestui semnal la o perioada ) Fie ()tx un semnal periodic, cu perioada T si (tx semnalul: 1 ⎧ TT ⎤⎡ () ,ttx ; −∈ ⎪ (10) ⎢ ⎥ () = 22tx ⎣ ⎦ ⎨ 1 ⎪ ,0 restin ⎩ Fie () ωX transformata Fourier a semnalului (tx si c coeficientii seriei Fourier 1)1k ) exponentiale asociati semnalului ( tx Atunci: (11) 21π () ,ωω , ZkkXc∈== k001 TT 6 Exemplu Se considera semnalul periodic din figura 1 () tx A t τ τ+T T−+− T 0Tτ = Fig 1 Semnal dreptunghiular, /1/4Tτ Coeficientii seriei Fourier exponentiale sunt: 3 τ (12) τω 2/sin1τωωjktjkk−−2/0 00 =⋅= eAdteAc⋅ k∫ 0 kT π iar puterea: 21(12`) τ⋅=AP X T Avand in vedere (7) rezulta: τ = AcA 00= T 2/sinkτ 0ω 22= = AcA (13) kk kπ ⎡ τωτω kkπ ⎞ ⎛ 00⎤ {} ,sinsgn1arg>⋅ 0 −+−==kcϕ ⎟ ⎜ kk ⎥⎢ 222⎠ ⎝ ⎦⎣ Tτ Spectrele de amplitudini si de faze sunt prezentate in figura 2 pentru cazul 4/1/= Fig 2a) Spectru de amplitudini pentru semnal dreptunghiular, /1/4Tτ = 4 Fig 2b) Spectru de faze pentru semnal dreptunghiular, /1/4Tτ= 7 Distorsiunile semnalelor Daca un semnal sinusoidal este aplicat la intrarea unui SLITC, atunci la iesire semnalul este sinusoidal In realitate, la iesire, semnalul este doar aproximativ sinusoidal, datorita neliniaritatilor inerente circuitelor reale Fie: ()() ψω+⋅=tAtxcos (14) 0 semnalul de la intrare si ∞ (15) () ( ) +⋅+=costkAAtyϕω ∑ 00kk 1k = semnalul de iesire Se defineste factorul de distorsiuni: 2 ∞ A∞ ⎛ k⎞ ⎜ ⎟ P∑ k∑ = ⎝ ⎠ 222kk= (16) d= Y= 112/PA Factorul de distorsiuni este radical din raportul dintre puterea continuta in armonicele superioare si puterea utila continuta in fundamentala semnalului 8 Desfasurarea lucrarii Se foloseste un generator de semnal, un osciloscop si un analizor de spectru Generatorul de semnal genereaza semnal periodic de forma sinusoidala, dreptunghiular, triunghiular Se aduce iesirea 50Ω a generatorului la intrarea osciloscopului (canalul 1) si la intrarea analizorului spectral 5 8 1 Se vizualizeaza si se deseneaza pe hartie milimetrica formele de unda ale unor semnale periodice generate de un generator de functii (dreptunghiular, sinusoidal, triunghiular) Semnalele sunt vizualizate in timp pe osciloscop si in domeniul frecventa folosind analizorul spectral Generator de semnal: frecventa fundamentala f0=1 MHz • amplitudine varf la varf App=1 V • pentru semnalul dreptunghiular: factorul de umplere 0 5 • Osciloscop: Scara de amplitudine 0 5 V/div; Scara de timp 0 5 μs/div • Analizor spectral: Center frequency: 5 MHz; Reference level: 10 dBm Span width: 1 MHz/div; RBW: 200 kHz 8 2 Sa se calculeze coeficientii Fourier pentru seria armonica Ak[V] in cazul ⋅ω (folositi identificarea), Semnalului sinusoidal () 00 5 sinxtV t= Semnalului dreptunghiular cu factor de umplere 0 5 (folositi exemplul din fig 1), Semnalului triunghiular Obs Semnalului triunghiular poate fi scris: 8111∞⎛⎞ ⋅ ω+ϕ= ω− ω+ ω− ω+ ()() ∑ 00 0000cos sin sin 3 sin 5 sin 7 kkxt A A k t t t t t=+ 2222⎜⎟ π ⎝⎠ 1357k= 8 3 Se masoara valorile efective ale armonicilor semnalelor de la punctul 8 1 cu un analizor de spectru Tab I Semnal sinusoidal F Uef [dBm] Uef[mV] U[mV] Ak[mV] f0 2f0 3f0 4f0 5f0 Tab II Semnal dreptunghiular, factor de umplere 0 5 Vom ignora armonicile de ordin par, ar trebui sa fie foarte mici (aproape de zero) f Uef [dBm] Uef[mV] U[mV] Ak[mV] f0 3f0 5f0 7f0 9f0 Tab III Semnal triunghiular f Uef [dBm] Uef[mV] U[mV] Ak[mV] f0 6 2f0 3f0 4f0 5f0 Se va folosi: VU [] ef 20 lg dBmU= [] ef refU Nivelul de referinta este 223 6mVU= ref UV [] VU= [] ef 2 8 4 Se reprezinta pe hartie milimetrica spectrele de amplitudini masurate si calculate Se compara valorile corespunzatoare 8 5 Calculati factorul de distorsiuni al semnalului sinusoidal considerand valorile efective din Tab I si relatia (16) 9 Prezentarea mediului Matlab MATLAB (MATrix LABoratory) este un program interactiv, dezvoltat de firma Math Works Inc , fiind destinat, în special, prelucrării numerice a datelor furnizate sub formă vectorială sau matriceală MATLAB integrează calculul numeric cu vizualizarea rezultatelor şi programarea, într- un mediu flexibil şi deschis dezvoltărilor ulterioare MATLAB este un sistem interactiv, al cărui element de bază este o matrice care nu necesită o dimensionare explicită Aceasta ajută în rezolvarea a numeroase probleme de calcul tehnic, în special pentru cele bazate pe formalismul matricial (şi vectorial), permiţând scrierea programelor într-un limbaj scalar MATLAB-ul a evoluat în ultimul timp, mai ales prin utilizarea sa în mediul universitar, ca instrument educaţional pentru cursuri de matematică, inginerie şi ştiinţă În egală măsură, MATLAB-ul se dovedeşte deosebit de eficient în cercetarea ştiinţifică şi dezvoltarea tehnologică în mediul industrial Moduri de lucru în MATLAB Lansarea în execuţie a programului MATLAB se face din WINDOWS După lansarea în execuţie, programul MATLAB intră în “modul de comandă”, afişând prompterul „ >> ” şi aşteptând introducerea unei comenzi de către utilizator Executarea unei comenzi este urmată, de obicei, fie de crearea unei variabile în spaţiul de lucru, fie de afişarea unui mesaj sau desenarea unui grafic De exemplu, comanda: >>v=0:10 7 va crea variabila v şi va afişa elementele acesteia pe ecran În afara modului de lucru “în linie de comandă”, în MATLAB se pot crea fişiere ce conţin instrucţiuni MATLAB, numite fişiere-M (deoarece au extensia “ m”) Un program MATLAB poate fi scris sub forma fişierelor “script” sau a fişierelor “function” Ambele tipuri de fişiere sunt scrise în format ASCII, iar algoritmul care a fost implementat poate fi urmărit cu foarte mare uşurinţă, dacă se cunosc convenţiile şi sintaxa MATLAB Aceste tipuri de fişiere, obligatoriu cu extensia “ m”, permit crearea unor funcţii noi care le pot completa pe cele deja existente Prin această facilitate, MATLAB-ul poate fi extins la aplicaţii specifice utilizatorului, care are posibilitatea să scrie noi proceduri Fişiere “script” Un fişier “script” este un fişier care conţine o secvenţă de comenzi MATLAB Prin apelarea numelului fişierului, se execută secvenţa MATLAB conţinută în acesta După execuţia completă a unui fişier script, variabilele cu care acesta a operat rămân în zona de memorie a aplicaţiei Aceste fişiere nu permit integrarea în programe mari, realizate pe principiul modularizării Fişierele script sunt folosite pentru rezolvarea unor probleme care cer comenzi succesive atât de lungi, încât ar putea deveni greoaie pentru lucrul în mod iterativ, adică în modul de lucru linie de comandă Fişiere “funcţie” Dacă prima linie a fişierului-M conţine cuvântul "function", fişierul respectiv este declarat ca fişier funcţie O funcţie diferă de un fişier "script" prin faptul că poate lucra cu argumente Variabilele definite şi manipulate în interiorul fişierului funcţie sunt localizate la nivelul acestuia Prin urmare, la terminarea execuţiei unei funcţii, în memoria calculatorului nu rămân decât variabilele de ieşire ale acesteia Fişierele funcţie sunt utilizate pentru extinderea MATLAB-ului, adică pentru crearea unor noi funcţii MATLAB Forma generală a primei linii a unui fişier funcţie este: function [param ieşire] = nume funcţie (param intrare) unde: function - este cuvânt cheie care declară fişierul ca fişier funcţie (obligatoriu); nume funcţie - numele funcţiei, adică numele sub care se salvează fişierul, fără extensie Nu poate fi identic cu cel al unui fişier-M preexistent param ieşire - parametri de ieşire trebuie separaţ i cu virgulă şi cuprinşi între paranteze drepte Dacă funcţia nu are parametri de ieşire, parantezele drepte şi semnul egal nu mai au sens param intrare - parametri de intrare trebuie separaţi cu virgulă şi cuprinşi între paranteze rotunde Dacă funcţia nu are parametri de intrare, parantezele rotunde nu mai au sens Semnale elementare 8 Impulsul unitate Impulsul unitate se defineşte în felul următor: 1, 0n= ⎧ nδ = [] ⎨ 0, 0n≠ ⎩ Utilizând proprietatea de deplasare în timp, putem scrie: ⎧ 01,nn= −= [] ⎨ 0nnδ 00,nn≠ ⎩ În MATLAB nu putem defini secvenţe de lungime infinită, trebuie precizat domeniul de valori pentru n Să se definească şi să se reprezinte grafic secvenţele: 1 [] = [] 1xnnδ =− 2 [] [ ] 21xn nδ 3 [] =+ pentru 10 10n−≤≤ [ ] 31xn nδ Toţi vectorii vor avea, deci 21 de elemente x1 % Generarea impulsului unitate clf; % sterge vechiul grafic n = -10:10; % generarea unui vector de la -10 la 10 d = [zeros(1,10) 1 zeros(1,10)]; % generarea impulsului stem(n,d); % reprezentarea grafică în timp discret xlabel('n');ylabel('Amplitudine'); title('Impulsul unitate'); axis([-10 10 0 1 2]); x2 n=-10:10; x2=zeros(size(n)); x2(12)=1; stem(n,x2),grid,title(‘x2[n]’),xlabel(‘n’) Momentul de timp n =1 corespunde celui de al 12-lea element al vectorului x3 n=-10:10; x3=zeros(size(n)); x3(10)=1; stem(n,x3),grid,title(‘x3[n]’),xlabel(‘n’) Momentul de timp n = -1 corespunde celui de al 10-lea element al vectorului Exerciţii Să se definească şi să se reprezinte graphic următoarele secvenţe: 9 =− pentru 120n≤≤ 1 [] [ ] 10 7 5xn nδ 2 [] = pentru 15 15n−≤≤ [] 10 6xnnδ Treapta unitate Treapta unitate se defineşte în felul următor: 1, 0n≥ ⎧ nσ = [] ⎨ 0, 0n< ⎩ Utilizând proprietatea de deplasare în timp, putem scrie: ⎧ 01,nn≥ −= [] ⎨ 0nnσ 00,nn< ⎩ În MATLAB nu putem defini secvenţe de lungime infinită, trebuie precizat domeniul de valori pentru n Să se definească şi să se reprezinte grafic secvenţele: 1 [] = pentru 10 10n−≤≤ [] 1xnnσ =− pentru 510n−≤ ≤ 2 [] [ ] 22xn nσ 3 [] =+ pentru 510n−≤ ≤ [ ] 32xn nσ x1 % Generarea treptei unitate n = -10:10; % generarea unui vector de la -10 la 10 u = [zeros(1,10) ones(1,11)]; % generarea treptei unitate stem(n,u); % reprezentarea grafica in timp discret xlabel('n');ylabel('Amplitudine'); title('Treapta unitate'); axis([-10 10 0 1 2]); x2 n=-5:10; x2=[zeros(1,7),ones(1,9)]; stem(n,x2),grid,title('x 2[n]'),xlabel('n') Momentul de timp n = 2 corespunde celui de al 8-lea element al vectorului x3 n=-5:10; x3=[zeros(1,3),ones(1,13)]; stem(n,x3),grid,title('x 3[n]'),xlabel('n') Momentul de timp n = −2 corespunde celui de al 4-lea element al vectorului Exerciţii Să se definească şi să se reprezinte grafic următoarele secvenţe: 1 [] = pentru 10 20n−≤≤ [] 10 7xnnσ 2 [] =− pentru 030n≤≤ [ ] 27xn nσ 3 [] =+ pentru 15 15n−≤≤ [ ] 31 8 3xn nσ 10